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1. Introduction 

The twisted mass formulation of Lattice QCD [1,2] has been studied extensively with Nf = 2 
dynamical flavours by the European Twisted Mass (ETM) collaboration. In this formulation of 
QCD, the Wilson term is chirally rotated within an isospin doublet. The effects of the strange 
and charm dynamical quarks are included through a mass-split doublet as discussed in [3, 4, 5, 6]. 
Results using two of the three lattice spacings discussed in these proceedings have recently been 
published in [7], where we describe our setup in more detail. Furthermore, at this conference, other 
investigations into the physics of A'y = 2 + 1 + 1 twisted mass fermions have also been presented: 
[8, 9, 10, 11, 12]. We will briefly describe our lattice setup and recapitulate our procedure for 
tuning to maximal twist and the strategy for the tuning of the heavy doublet in section 2. We give 
an overview of the runs we have earned out in section 3, where we also examine the status of the 
tuning. Finally section 4 gives preliminary results for some observables in the light-quark sector, 
obtained using fits to next-to-leading order (NLO) SU (2) chiral perturbation theory. 

2. Lattice setup 

In the gauge sector we use the Iwasaki gauge action [13] since it improves the behavior of 
the lattice theory in relation to the unphysical first order phase transition for values of the hopping 
parameter K around its critical value K^nt (see [7] and references therein). With this gauge action 
we observe indeed a smooth dependence of phase sensitive quantities for K ~ JQnt. The fermionic 
action for the light doublet is given by: 

Si = a^Y. i^/ (^) [^] + '"0./ + T3] Xi (x) } , (2. 1) 

using the notation used in [7]. In the heavy sector, the action becomes: 

Sh = a'^Y. (-^) [Dw[U]+ mo,h + il^a 75 + Ms 'l^s] Xh {x) } . (2.2) 

X 

At maximal twist, physical observables are automatically d'{a) improved without the need to de- 
termine any action or operator specific improvement coefficients. The gauge configurations are 
generated with a (Polynomial) Hybrid Monte Carlo updating algorithm [14, 15, 16], where the 
HMC is used for the light doublet and the PHMC for the heavy doublet. 

Tuning to maximal twist requires to set mo,/ and mo,/i equal to some proper estimate of the crit- 
ical mass merit = "icrit(j8) [3]. As has been shown in [4], this is consistent with i^{a) improvement 
defined by the maximal twist condition ampcAC.i = (see also ref. [7]). The numerical precision 
at which the condition mpcAC./ = is fulfilled in order to avoid residual large tff{a^) effects when 
the pion mass is decreased is, for the present range of lattice spacings, |£//i/| ^0.1, where e is 
the deviation of mpcAC,/ from zero [17, 18]. As explained in [7], tuning to K^rit was performed 
independently for each set of values of /X;, jXa and Us- From table 1 we observe that the estimate 
of K"crit depends weakly on /i/. The heavy doublet mass parameters jtttj and jJ-g should be adjusted 
in order to reproduce the values of the renormalized s and c quark masses. The latter are related to 
Ha and via [3]: 

1 Zp 

{nh,c)R = —il^oT—l^s), (2.3) 



2 



Light hadrons from Nf = 2 + 1 + 1 dynamical twisted mass fermions 



S. Reker 



where the minus sign corresponds to the strange and the plus sign to the charm. In practice we fix 
the values and Hg by requiring the resulting kaon and D meson masses to match their physical 
values. A detailed description of the determination of the kaon and D meson masses has recently 
been given in [19]. 



3. Ensemble overview 

We list in table 1 the action parameters for the runs considered in our current analysis. Those 
runs labeled with an asterisk (*) ai^e ongoing at the time of the writing of these proceedings and 
have incomplete statistics, all other runs have around 5000 thermalized trajectories with length 
T = 1 . Ensemble names which end in s or c are used to control the tuning of the strange and charm 
quark masses respectively. This is not an exhaustive overview of all runs performed within of our 
Nf = 2 + 1 + 1 work. Other runs have been canied out in the context of tuning, reweighting, finite 
volume effects analysis and in order to measure the renormalization factors. 



Ensemble 


/3 


^cht 


<3jU, 






{L/a)^ xT/a 


A30.32 


1.90 


0.1632720 


0.0030 


0.150 


0.190 


32^ X 64 


A40.32 




0.1632700 


0.0040 






32^ X 64 


A50.32 




0.1632670 


0.0050 






32^ X 64 


A60.24 




0.1632650 


0.0060 






24^ X 48 


A80.24 




0.1632600 


0.0080 






24^ X 48 


A100.24 




0.1632550 


0.0100 






24^ X 48 


A80.24S 




0.1631204 


0.0080 


0.150 


0.197 


243 X 48 


A100.24S 




0.1631960 


0.0100 






243 X 48 


B25.32 


1.95 


0.1612420 


0.0025 


0.135 


0.170 


32^ X 64 


*B35.48 




0.1612400 


0.0035 






48^ x96 


B35.32 




0.1612400 


0.0035 






32^ X 64 


B55.32 




0.1612360 


0.0055 






32^ X 64 


B75.32 




0.1612320 


0.0075 






32^ X 64 


B85.24 




0.1612312 


0.0085 






24^ X 48 


*D1 15.64 


2.10 


0.1563640 


0.00115 


0.120 


0.1385 


64^ X 128 


D15.48 




0.1563610 


0.0015 






48^ x96 


D20.48 




0.1563570 


0.0020 






48^ x96 


D30.48 




0.1563550 


0.0030 






48^ x96 


*D45.32sc 




0.1563550 


0.0030 


0.0937 


0.1077 


32^ X 64 



Table 1: Summary of the A^f = 2 + 1 + 1 ensembles generated by ETMC at three values of the lattice 
coupling j3 = 1.90, j3 = 1.95 and j3 =2.10. From left to right, we quote the ensemble name, the value of 
inverse coupling j3, the estimate of the critical value Kent, the light twisted mass a/i/, the heavy doublet mass 
parameters ajXc, and ajig and the volume in units of the lattice spacing. Our notation for the ensemble names 
corresponds to X/i/.L, with X referring to the value of j3 used. 
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3.1 Tuning to maximal twist 

Figure 1 shows the status of the tuning for the main ensembles considered in these proceedings. 
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Figure 1: Status of the tuning. The ratio wjpcac,;/M/ is plotted as a function of the mass parameter 2Z?o/i/. 
When |topcac,//M/| ^0.1, the ensemble is adequately tuned. Orange, blue and green symbols respectively 
correspond to j3 = 1.90, j3 = 1.95 and /3 =2.10 ensembles respectively. 



3.2 Heavy doublet tuning 

Figure 2 shows the dependence of {2m^ — mp^) and mj) upon the light pseudoscalar mass 
squared for both ensembles, together with the physical point denoted by the black star. The kaon 
mass appears to be properly tuned at j3 = 1.95. The ensembles at j8 = 1.90, a^g = 0.190 and 
/3 =2.10 appear to have a value of the strange quark mass lai^ger than the physical one, while the 
red point at /3 = 1.90, ajJ-s = 0.197 appears to be well tuned. The D meson appears heavier than 
in experiment for all three values of the lattice spacing. We currently have runs ongoing at both 
j3 = 1.90 and j8 =2.10 with lower charm masses. 

4. Light meson chiral perturbation theory fits 

In order to extract the lattice spacing and light quark mass from our data, we perform a NLO 
SU (2) chiral perturbation theory fit of the m,: and fz lattice data. We group our ensembles into sets 
with the same lattice spacing (set A at j3 = 1.90, S at jS = 1.95 and set D at j8 = 2.10). We have 
performed fits for various combinations of these sets, using the procedure described in [7]. We use 
continuum formulae and currently correct for finite volume effects as described in [20]. Note that 
since the quark mass enters the ^PT expression, in order to combine ensembles at different lattice 
spacings, we need to know the renormalization factor of the quark mass = 1 /Zp, a computation 
which is not yet complete. Assuming that Zp is effectively a function of /3 in the range of parameters 
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Figure 2: 2ot|. — mpg and mo as a function of mpg . The physical point is shown (black star) [21]. Data points 
have been scaled with the lattice spacing a = 0.0863(4) fm for j3 = 1 .90, a = 0.0779(4) fm for ^ = 1 .95 and 
a = 0.0607(2) fm for j3 = 2.10, where the errors quoted on the lattice spacing are only statistical. 



we are considering, we can fit the ratio of those Zp-values and lattice spacings and extract lattice 
spacings from the combined fit. In every fit we use as inputs the physical values of fji and m„, and 
extract /o, h, U and the lattice spacing. The results are listed in table 2, while figure 3 shows the fit 
to sets A, B and D combined. 



set 


pts 


/o(MeV) 


I3 


h 


«/3=i.9o(fm) 


«i3 = 1.95(fm) 


a/3=2.io(fm) 


A 


5 


120.96(7) 


3.44(6) 


4.77(2) 


0.0859(5) 






B 


4 


121.15(8) 


3.70(7) 


4.67(3) 




0.0782(6) 




A&B 


9 


121.03(5) 


3.54(5) 


4.74(2) 


0.0861(4) 


0.0778(4) 




A&D 


7 


120.99(7) 


3.42(7) 


4.76(3) 


0.0861(5) 




0.0606(3) 


B&D 


6 


121.20(8) 


3.68(7) 


4.65(3) 




0.0785(6) 


0.0609(3) 


A&B&D 


11 


121.05(5) 


3.53(5) 


4.73(2) 


0.0863(4) 


0.0779(4) 


0.0607(2) 



Table 2: Results from the NLO SU{2) ;ifPT fits for various combinations of the ensembles. Errors are 
statistical only, extracted from 200 bootstrap samples. The column "pts" refers to the number of ensembles 
used in that fit. 



The observed agreement between the extracted parameters suggests that our data for mps and 
fps are fairly well described by NLO SU (2) chiral perturbation theory. Using the spread of pa- 
rameters as a rough estimate of the systematic error, it appears to be smaller than the statistical 
error for all quantities, with the exception of Z3. A more complete analysis of the systematic effects 
(analogous to [22]) is in progress. We expect to extend our analysis by including twisted mass chi- 
ral perturbation theory formulae as described in [23] and to use twisted mass finite volume effects 
formulae [24] when our neutral pion measurements are more complete. 
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Figure 3: The charged pseudoscalar mass ratio m^^/2Bo^i and the pseudoscalar decay constant /ps as a 
function of the mass parameter IBofii, for the combined ensembles at j3 = 1.90, )3 = 1.95 and j3 = 2.10. 
The scale is set by a/iphys, the value of ajj.) at which the ratio /ps °°V'"£ assumes its physical value [21] 
fn/niji = 130.4(2)/ 135.0 (black star). Open symbols refer to runs with full statistics, but not properly tuned 
to maximal twist within our criterion. Runs not at full statistics and those aimed at controlling the tuning of 
the strange and charm mass are not included in the plot. 



5. Summary and Outlook 

We have given an update of the status of the runs performed by the ETM Collaboration using 
Nf = 2 + 1 + 1 flavours of Wilson twisted mass fermions. We have given first results at a new finer 
(jS = 2. 10) lattice spacing and attempted to combine them with existing datasets at two other lattice 
spacings (jS = 1.90 and j8 = 1.95). The production of ensembles at the finest lattice spacing is still 
ongoing. As already stated, a complete control of the different systematic effects present in chiral 
fits of pion observables is still missing. ETMC is currently pursuing the direct determination of the 
renormalisation factor Zp [12] appearing in the fits combining ensembles at several values of the 
lattice spacing. The preliminary results presented in this work are nevertheless very encouraging 
and suggest a fairly good description of our lattice data for mps and /ps by NLO SU (2) chiral 
perturbation theory. 
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